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We reconsider some general aspects about the mean field thermodynamical description of the 
astrophysical systems based on the microcanonical ensemble. Starting from these basis, we devote a 
special attention to the analysis of the scaling laws of the thermodynamical variables and potentials 
in the thermodynamic limit. Geometrical considerations motivate a way by means of which could be 
carried out a well-defined generalized canonical-like description for this kind of systems, even being 
nonextensive. This interesting possibility allows us to extend the applicability of the Standard 
Thermodynamic methods, even in the cases in which the system exhibits a negative specific heat. 
As example of application, we reconsider the classical Antonov problem of the isothermal spheres. 
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I. INTRODUCTION 



II. MICROCANONICAL DESCRIPTION 



Traditional thermodynamics is not able to describe 
the astrophysical systems because of their nonextensive 
nature: They are nonhomogeneous systems whose total 
energy does not scale asymptotically with the particles 
number N during the increasing of the system size, as 
well as they exhibit energetic regions with a negative 
specific heat, which makes the canonical description un- 
suitable for describing their thermodynamical properties 

QiiiiiHiillS. 

An adequate thermodynamical description for astro- 
physical systems can be carried out starting from micro- 
canonical basis under certain conditions l^, This 
possibility is very attractive because the microcanoni- 
cal ensemble is the only well-defined statistical ensem- 
ble: since dynamics is always well defined, and because 
microcanonical ensemble is a "dynamical" ensemble. 

The aim of the present paper is to reconsider some 
general questions about the thermodynamic formalism 
for astrophysical systems based on the consideration of 
the microcanonical ensemble. The organization of the 
paper is the following. Sec. II is devoted to review some 
aspects of the microcanonical mean field description of 
the astrophysical systems in order to provide a general 
background for the analysis of the scaling laws of the 
thermodynamical variables and potentials. Some details 
about the Legendre formalism are also presented. Impli- 
cations of the geometrical aspects of the microcanonical 
ensemble on the applicability of a generalized canonical- 
like description are discussed in Sec. III. Finally, in Sec. 
IV some conclusions are drawn. 



A direct microcanonical calculation in a selfgravitating 
system is an formidable task which can be performed by 
using Montecarlo numerical simulations (see for example 
in refs.0, ^3)- The statistical description of A^-body 
selfgravitating system is usually carried out by using a 
mean field approximation when the particles number N 
is large enough |^ . We will perform in this section 
a brief review about some aspects of this microcanonical 
mean field formalism for the astrophysical systems. 



A. The microcanonical mean field approximation 

Let us to consider an astrophysical Hamiltonian sys- 
tem composed by a huge number N of identical particles 
enclosed in a rigid container with a characteristic linear 
dimension L. The consideration of this rigid container 
is a regularization procedure for the thermodynamical 
description for this kind of system which avoids the par- 
ticles evaporation 0. Let us also suppose that these 
particles interact among them by means of the gravity 
and short-range interactions, as example, forces with an 
electrostatic origin. The Hamiltonian of this system can 
be discomposed into two terms as follows: 



H 



N 



In ■ 



V, 



(1) 



where Vn is the potential energy associated to the New- 
tonian gravitational interaction, while /i^r contains the 
energy contribution of the short-range interactions and 
the kinetic energy of the particles. 

The entropy function for the selfgravitating system 
in the microcanonical description is obtained from the 
Boltzmann principle: 
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where the microcanonical accessible volume W {E, N; L) 
is given by: 



W {E, N; L) = Sp [5 {E ~ Hn)] Seq, 



(3) 



being 6eQ a suitable energy constant which makes W 
dimensionless. Quantum effects or/and the consider- 
ation of a natural size for the particles will be taken 
into account in an implicit manner in the equation 
which is a necessary consideration in order to avoid the 
short-range singularity of the Newtonian interaction, and 
therefore, the gravitational collapse 0, . 

We are interested in describing the large N limit. This 
aim can be carried out by using the following proce- 
dure: we partition the physical space in cells {cq,} whose 
characteristic linear dimensions are much larger than the 
effective radio of the short-range interactions, but very 
short in comparison with the characteristic linear dimen- 
sion L of the container. This assumption supposes that 
this container is very large. We denoted hy Ua = n (ja) 
the number of particles inside the cell Cq., being the 
position of its center. 

The short-range interactions among the particles be- 
longing to different cells can be neglected, and therefore, 
the term hpj of the Hamiltonian Q can be approximated 
by: 



(4) 



where /i„^ is the energetic contribution of the Ua parti- 
cles inside the cell Cq. Hereafter, we refer this energetic 
contribution as the internal energy. The potential en- 
ergy Vn associated to the gravitational interaction can 
be approximated as follows: 



Vr 



N 



V{n] = -Y: 



Q>/3 



fa - r^l 



(5) 



Taking into consideration all approximations exposed 
above, the microcanonical volume W (jSJ can be rewritten 
as follows: 



W{E,N,L)^ J2dD iN~Y.na]w[n;E,L], (6) 



{no} 



being: 



The functional W [n; E, L] is given by: 



W[n;E,L] = 




<5eo<^ \E~V[n]~ ^e^. 



X exp 



(8) 



where Sa {ea,najVa) — In [cj^ (cq., na, is the Boltz- 
mann entropy associated to the subsystem of particles 
inside the cell Cq whose internal energy is Bq, being Va the 
volume of this cell. The quantity cJq, (ec, Wq) is the 
number of microstates of this subsystem with the above 
macroscopic configuration: 



(9) 



where Sea is a suitable energy constant. Since /i„^ in- 
volves only short-range interactions, this subsystem will 
look-like as an extensive system for Ua large, and there- 
fore, the entropy per volume Sq depends only on the in- 
ternal energy density = ea/^a and the particles den- 
sity Pa — ria/va as foUows: 



Sa {CcK 7 ^ct ; '^a ) /'^a — ^ {^a i Pa ) 



(10) 



The microcanonical volume W can be conveniently 
rewritten by using the following mean field approxima- 
tion: 



W ^ Wmf {E, N,L) = aJvp (r) Ve (r) exp {S [e, p]) x 

xS{N-N[p])S{E^H[e,p]), (11) 

where A is an unimportant constant which involves the 
constants Se^s and VaS. The functionals 



S[e,p]^ / dhs{eiv),piv)}, 



H[e,p, 



(12) 



d^r e (r) + -mp (r) (r) , (13) 
z 



N[p]^ f d'rpir), (14) 

represent the total entropy, energy and particles number 
for a given profile with e (r) and p (r) , being (j) (r) the 
Newtonian potential: 



I (r) = g [p; r] 



Gmp (ri) d^Ti 
|r-ri| 



(15) 
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where Q [p; r] is the Green solution of the Poisson prob- 
lem: 



thermodynamical variables of the system are scaled by 
some scaling parameter a as follows: 



A(j> = AnGmp. 



(16) 



The expression of the microcanonical volume in mean 
field approximation can be rewritten in order to avoid 
the p dependence of the Newtonian potential as follows: 

Wmf {E, N,L)=aJvp (r) Ve (r) Vcj) (r) exp (5 [e, p]) x 



N - 
E - 
V ■ 



N{a) 
E{a) : 
V{a) 




W{a) =exp {a\nW). 

,(21) 

In analogy with the extensive properties of the traditional 
systems, we will analyze the necessary conditions for the 
existence of the following power law self- similarity scaling 
behavior of the microcanonical variables E, N and L for 
the astrophysical systems: 



x5{cl>iv)~-g[p;v]}SiN-N[p])SiE-H[e,p, <!>]). 

(17) 

Wmf{E,N,L) can be rewritten again by using the 
Fourier representation of the delta functions, yielding: 



Wmf iE,N) 



oo p-\-oo 



dkd7] 
(2^ 



VpVeDclfDhx 



exp {£ [e,p,(/);z,zi, J]}, 



(18) 



where z = (3 + ik and zi — p -\- irj with /5, ?/ e R, being 
the functional C [e, p, 0; z, zi, J] defined by: 



C [e, p, 0; z, zi, J] = 5 [e, p] + z{E ~ H [e, p, </>]) + 
+zi (7V-iV[p]) + J.(0-^[p]). 



(19) 



The functional term J ■ {4>~Q[p]) — 
/g3 d^T J (r) {0 (r) — Q [p; r]} appears as consequence 
of the Fourier representation of the delta functional 
5 {cf){r) — Q [p]r\}. Here, J (r) is a complex function, 
J (r) i (r) + i/i (r), with j (r) e M. 



B. Scaling properties of the thermodynamical 
variables and potentials 

The integration of expression (|18(l is usually carried out 
by using the steepest decent method. The application of 
this method is based on the asymptotic behavior of the 
thermodynamical variables and potentials in the many 
particle limit ^ 1. 

The presence of an additive kinetic part in the Hamil- 
tonian of certain system leads to an exponential growing 
of the microcanonical volume W with the N increasing, 
and therefore, the Boltzmann entropy will grow propor- 
tional to N in the many particles limit, 5*3 = In cx A^. 

The usual thermodynamic limit for the extensive sys- 
tems: 



E N 

N ^ oo, keeping constant — and — , (20) 

where V is the volume of the system, is directly related 
with the extensive properties of these systems when the 



N 

E - 



E{a)-- 



L{a) = a^^L, 



W^W{a)^ exp(alnTy), 
(22) 



where A^; and are certain constant scaling exponent 
which lead to an extensive character of the Boltzmann 
entropy. This kind of self-similarity behavior is directly 
related with a thermodynamic limit of the form: 



iV — > cx), keeping constant 



E 



and 



(23) 



The existence of this kind of self-similarity condition 
allows a considerable simplification of the thermodynam- 
ical description: the study can be performed by setting 
A'^ = 1 and considering the A'^-dependence in the scaling 
laws by taking a — N. This scaling behavior is very use- 
ful in numerical experiments, since it allows us to extend 
the results of this kind of study on a finite system to much 
bigger systems. Contrary, the nontrivial A^-dependent be- 
havior of the thermodynamical variables and potentials 
leads to a complication of the analysis. 

In order to satisfy this scaling behavior for the global 
variables E, N, L and the Boltzmann entropy 5*^, the 
local functions e(r), p{r), (f>{r) and s{e,p;(l)) should be 
scaled as follows: 



p^ p{a)=a^''p} 



0(a) 
s (a) 



(24) 



Since the characteristic particles density behaves as 
Pc ^ N/L^, the scaling exponent for the particles density 
is Ap = 1 — 3Ai. From the expression of the Newtonian 
potential H15|l is derived that its characteristic unit is 
(f)c ^ PcL^, and therefore, A^ = 1 — A^. The energy scal- 
ing exponent is equal to the scaling exponent of the total 
gravitational potential energy, so that, A^; = 2 — A/,. The 
other scaling exponents are obtained by using identical 
reasonings. All these scaling exponents depend on the 
scaling exponent A^ as follows: 



Ap = 1 - 3Al = As, A^ = 1 - Al, 
Ae = 2-4Al, Ae = 2- Al. 



(25) 
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In order to satisfy these scaling laws is also necessary 
that the entropy density exhibits to the following scaling 
behavior: 



s{a^'e,a^''p) ^a^<'s{e,p), (26) 

This scaling property is satisfy if s (e, p) obeys to the 
following functional form: 



s(e,p) = pF(e/p") 



(27) 



where rj — Ag/Ap, being F (x) an arbitrary function. It 
is easy to show that the functional form l|27|l leads to the 
following relation between the pressure p and the internal 
energy density e: 



P = 7e, 



(28) 



where 7 = 77 — 1. There are some well-known Hamilto- 
nian systems which satisfy this kind of relation, as ex- 
ample, the system of nonrelativistic or ultrarelativistic 
noninteracting particles, without matter if they obey to 
the Boltzmann, Fermi-Dirac or Bose-Einstein Statistics. 
The scaling parameter A^ and A^ are obtained from the 
parameter 7 as follows: 



the functional integral (|18|l will come from the maxima 
of the functional £ [e, p, 0; z, zi, J]. This last one can 
be rephrased as the Legendre functional, in which the 
entropy functional S [e, p] is maximized under the con- 
strains of the energy and particles number, where it is 
also taken into account the relation between the New- 
tonian potential and the particles density. As elsewhere 
discussed lOj, the existence of a multimodal functional 
C [e, p, (j); z, zi, J] for a given value of the total energy E 
tells us about the existence of several quasi-equilibrium 
profiles, which could be related with the existence of 
phase transitions [loL ll^ . When there is only one sharp 
peak, the Boltzmann entropy can be appropriately esti- 
mated as follows: 



Sb {E, N, L) 



max 

e, P, 4> 



min C[e,p,(j);P,p,j] 



(31) 



where P and p are the canonical parameters (real num- 
bers), while j is a real field which acts as a Lagrange 
multiplier for the relation between the Newtonian poten- 
tial and the particles density. The stationary conditions 
lead to the following relations: 



13 = d^s (e,p) , p 



-^Pm(t)-g [j]+dpS (e, p), j = ^l3mp, 

(32) 



At 



1 



37 - r 



57-1 
87-1' 



(29) 



This result evidences that the power laws form for the 
self-similarity conditions (|22l) can be only satisfied by a 
reduced group of models whose microscopic picture obeys 
to the relation H28|) . Since 7 = | for the ideal gas of 
particles, the scaling exponents for the Antonov problem 
Q and the selfgravitating fermions model [isl IT^ are 



given by A_e = | and A^ 



and therefore, they obey 



to the following thermodynamic limit: 



N ^ 00, keeping constant 



E 



and LN^ . 



(30) 



This thermodynamic limit was established in the ref.|l(]| 
for the self-gravitating nonrelativistic fermions by using 
other reasonings. Since the selfgravitating relativistic gas 
and classic hard sphere model J^J do not obey to the 
relation (|28|l . they do not posses a scaling behavior with 
a power law form (|22|l . 



C. Legendre formalism 

For the sake of simplicity, let us consider firstly those 
models exhibiting a power law thermodynamic limit (|23|l . 
In such cases, the steepest decent method leads to a 
thermodynamic formalism very analogue to the one used 
for the extensive systems. The main contribution of 



N[p]^N, H[e,p,<l>]^E, <l> = g[p]. 



(33) 



The first relation of (|32|l states that the canonical pa- 
rameter (3 is constant through all points of the system. 
Therefore, it is convenient to use the Planck density 
Pc iP, p) = max [(3e — s (e, p)] instead of the entropy den- 

€ 

sity, which allows us to rewrite the relations H32|l as fol- 
lows: 

e = dfjpc (/?, p), p^ -fSmcj) ~ dpPc (/3, p) , (34) 
where the field j was eliminated and the identity 



dpS (e, p) = -dpPc (/?, p) 



(35) 



was taken into account. The second relation in the equa- 
tion H34|) allows the obtaining of the state equation for 
the particles density, p = p (0; /3, p). The consideration of 
the constrains 1)33(1 lead to a self-consistent integral equa- 
tions system whose solution represents the most probable 
equilibrium configuration of the astrophysical system for 
a given total energy E. 

Fluctuations around this equilibrium profile for the 
particles density p could be estimated by using a Gaus- 
sian approximation as follows: 



(36) 



eq 



5 



Since £ grows proportional to N when A'' ^ 1, the rel- 
ative fluctuations depends on the system size as Sp/p ~ 
1 / VN, in analogue form that the traditional systems do 
obey. Although all reasonings have been made for as- 
trophysical models with a power law thermodynamical 
limit (|23|l . the estimation of the relative fluctuations of 
the particles density p (|36|l evidences the general appli- 
cability of Legendre formalism exposed above. This con- 
clusion is straightforward followed from the vanishing of 
the relative fluctuations of the particles density p due to 
the general linear growing of the functional C with the 
N increasing. 

This theoretical fact explains why several tools of the 
Traditional Thermodynamics can be extended for the as- 
trophysical objects in spite of their nonextensive charac- 
ter. This conclusion has been elsewhere experimentally 
confirmed by observations. For example, the evidence 
about the isothermal character of the core in globular 
clusters and elliptical galaxies, where these last ones dis- 
play a quasi-universal luminosity profile described by de 
Vaucouleur's i?^/^ law 0|; the presence of a dark mat- 
ter halo density profile decreasing as at large dis- 
tances in the spiral galaxies, which is is also related with 
the isothermal distributions However, in spite of 

the common applicability of several tools of the standard 
thermodynamic formalism, the thermodynamical prop- 
erties of astrophysical systems turn to be very different 
from the ones exhibited by the ordinary extensive sys- 
tems due to the long-range character of the gravitational 
interaction, which leads to the existence of some striking 
phenomena, i.e. the gravitational collapse and the ex- 
istence of a negative heat capacity for certain energetic 
region Q, i S |l i i .S, Ofli] 



III. GEOMETRICAL CONSIDERATIONS 

Is it only microcanonically that could be performed 
a well-defined description for the astrophysical systems? 
or, will there be also a well-defined canonical-like descrip- 
tion? In our opinion, the answer of the last question is 
yes, which will be proved in the following subsections. 



A. Physics in the microcanonical ensemble is 
reparametrization invariant. 

Let us consider certain functional (f = (f {E, N) which 
is a bijective function of the total energy E (there is a bi- 
jective map between E and ip) . In this case, the function 
(fi {E, N) is also an integral of motion for the system. Let 
us introduce the microcanonical state density fi^: 



((^ {E, N),N)= {E, N)-^ [Hn {X) , N]] dX. 



By using the identities: 

S[^{E,N)-ipiHN,N)] = 
and 



dip iE,N) 



dE 



S [E - Hn) , 
(38) 



dip iE,N) 



dE 



rtE (E, N) , (39) 



is straightforward derived the reparametrization invari- 
ance of the microcanonical probabilistic distribution func- 
tion: 



5[p~P>{Hn,N)] = 



1 



{E, N) 



5{E-Hn) 



LOM {X; p, N) = UJM {X- E, N) 



(40) 



From this property immediately follows the 
reparametrization invariance of the microcanonical 
expectation values: 



F (X) UJM {X; p,N)dX= F (X) ujm [X] E, N) dX. 



Fm{^,N)^Fm {E,N). 



(41) 



The previous results point out the following conclu- 
sion: the microcanonical description is equivalently per- 
formed by using the thermodynamical variables [E, N) 
or {p, N). We say that {E, N) and {p, N) are two equiv- 
alent representations for the abstract space of the mi- 
crocanonical macroscopic description ^n. The set of 
all those representation changes among equivalent rep- 
resentations of Sat constitute the group of diffeomor- 
fisms or reparametrizations of Qat, which is denoted by 
Diff (3jv) . Thus, Physics in microcanonical ensemble is 
invariant under the reparametrizations changes o/ Sat. 



B. Extensivity of the Planck thermodynamic 
potential. 

The analysis of the necessary conditions for the en- 
semble equivalence starts from the consideration of the 
Laplace transformation between the microcanonical and 
the canonical-like partition functions, fl^ {p, N) and 



(37) 



{f3^, N) = / exp {-(3p) {ip, N) dp. 



(42) 
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Let us to show that it is always possible to choose a 
representation {ip,N) of S^v where the Planck potential: 



(43) 



is extensive for a system with a power law self-similarity 
scaling behavior (|22|l . 

Let W be the microcanonical phase-space accessible 
volume: 



W (E, N)^n (E, N) Se, 



(44) 



For continuous variables, W is only well-defined after a 
coarsed grained partition of phase space, which is the 
reason why it is considered a small energy Se in order to 
make W dimensionless. However, during the thermody- 
namic limit TV — !■ cx), we can make an arbitrary selection 
of Se, whenever this last be small. In this case, the micro- 
canonical phase-space accessible volume W appears as a 
reparametrization invariant function: 



The simplest choice for (p{E,N), in accordance with 
the thermodynamic limit (|23|l . is given by the scaled en- 
ergy £ : 



(49) 



In general, the extensivity of ip {E, N) is ensured by con- 
sidering the following dependence: 



ip{E,N) = N(t> {£/N) , 



(50) 



where (j) (e) is a bijective function of e. All those rep- 
resentation changes </? G Dif f {'^m) preserving the ex- 
tensivity of £ in the thermodynamic limit constitute a 
subgroup, which can be referred as the homogeneous sub- 
group Mc C Diff (Qat). All the admissible representa- 
tions for the canonical-like description are parametrized 
by considering all transformations of A^c on the scaled 
energy £, the pair {£,Mc)- 



W = nsSe = n^S(p, 



(45) 



where Se and Sip are very small in order to satisfy the 
condition: 



Se 



dip 



BE 



Sip. 



In this case, the Boltzmann entropy: 

Sb = InVF, 
becomes a scalar function: 

Sb [p, N) = Sb {E, N) . 



(46) 



(47) 



Taking into account all the exposed above, the equa- 
tion can be rewritten as follows: 



exp [-V^ (/3^, 



N)] = j exp [-/3^^ + Sb (p, N)] ^. 



(48) 

Sb {ip, N) is extensive when the thermodynamic limit 
(|23|l is carried out. The extensivity of the Planck po- 
tential Vip (Pip, N) when N is tended to infinity is guar- 
anteed by choosing a functional p {E, N) which must be 
extensive when the thermodynamic limit is performed . 
This is a necessary condition for the ensemble equiva- 
lence. Although in the microcanonical description any 
representation {p>,N) of Sat can be considered as ther- 
modynamic variables, in the canonical-like description 
are only admissible those representations satisfying the 
extensity condition. This condition is arisen as the gen- 
eralization of the additivity condition of the traditional 
Thermodynamics. 



C. Ensemble equivalence 



When N is tended to infinity by using a representation 
TZ^ e {£,Mc) in the equation the main contribu- 
tion of this integral will come from the maxima of the 
exponential argument. The equivalence between the mi- 
crocanonical and the generalized canonical ensemble will 
take place when there is only one sharp peak. In such 
cases the validity of the Legendre transformation is en- 
sured: 



V^{l3^,N)^Tmn[l3^ip- Sb{p},N)]. (51) 
The minimization leads to the conditions: 



d_ 

dp> 



Sb{v,N),K^ 



dpi- 



■5s((^,7V) < 0. (52) 



In this case, = d^fi^p appears as a generalization 
of the curvature tensor of the microcanonical thermosta- 
tistical theory of Gross [l3| . The topology of this tensor 
defines the ordering information of the system by using 
a microcanonical description. As usual, all those points 
of Sat where > can not be accessed by using the 
canonical-like description in the ip-representation, TZp. 
However, it can be easily shown that the sign of is 
representation dependent. 

Let V™ {p, N) be the microcanonical Planck potential: 



v::{p,N) 



^g^Sip,N) 



Sb (<^, N) . 



(53) 



V™ [p, N) becomes in the canonical Planck potential 
when the ensemble equivalence takes place in the thermo- 
dynamic limit. Let us now consider two representations 
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TZip^ and TZip^ of (5,A^c)i where the reparametrization 
change Ti^2 ^20^1^ G A^c '■ 

ip2^ {V20(Pi^){ipi) ^ 'P2{Vl,N) . (54) 

Since preserves the extensivity of Lpi, hence, 

Lp2 (vi, N) is a homogeneous function: 



as well the A'^-independent thermodynamic variables e 
and ip = ^/N. In this case, the canonical parameter /3 
and the Af-independent curvature k are denoted by: 

= d,s, h = (62) 
in the TZg representation, or 



and therefore, it obeys to the identity: 



dip2 , ,^dip2 



(55) 



(56) 



Taking into consideration the definitions (|52|) and (|53|l . 
as well as the identity (|56|l . it is straightforward derived 
the identities: 



ckpi 
dip2 



difi 
dip2 



V"\{^2,N)^V'^Avi,N)+N 



d(p2 
dN 



(57) 



(58) 



(59) 



It is evident from H58|) that the sign of is 



non invari- 
ant under the reparametrization changes A4c- Note that 
(3^ obeys the transformation rule for covariant vectors in 
a Riemannian geometry. However, K^p does not obey 
the correct transformation rule of second-rank covariant 
tensors. This is the reason why the curvature tensor is 
not a real tensor, and the ordering information which it 
contains is non invariant to under the reparametrization 
changes of A^c- A question arises: is it possible to choose 
a representation TZ^ in which the ensemble equivalence 
takes place? The answer is yes, and we will prove it now. 

As already said, when it is used the scaled energy rep- 
resentation TZs = {£, N) H49|) . there is an energetic region 
where the astrophysical systems exhibits a negative heat 
capacity: Kg > 0. Let us consider another representation 
TZq, = (<&, TV) which can be parametrized by a bijective 
function (p {e) as follows: 



$ = Nip{e) with e = £/N. 



(60) 



In order to disregard the A'^-dependence, let us work with 
the Boltzmann entropy per particle in the thermody- 
namic limit: 



s (e) — lim 



SBi£,N) 
N 



(61) 



(63) 



in the 7?.$ representation. The notation d^s is equivalent 
to ds/de. The expression H58|l can be rewritten in this 
case as follows: 



dpPp (^) = {d,^ {e))-' [d,P, (6) - /?, (e) a. In (9,^ (e))] 

5^/3^ = -(5,^(e))-'a(e), (64) 
where the function a (e): 

d,P,-MMid,ipie)) = -a{e). (65) 

should be positive in order to ensure the ensemble equiv- 
alence in the 7?.$ representation. Direct integration of 
(|^ yields: 



de(p (e) = CPe (e) exp 



de 



(66) 



where C is a positive constant which could be set as unity. 
It is easy to see that a convenient choice for the function 
a (e), from the theoretical viewpoint, could be given by: 



-k, (e) if k, (e) < 0, 

t (e) = <' /?2 (e) = /32 ^ const if k, (e) = 0, 
fee (e) if k, (e) > 0, 



(67) 



which ensures the positivity of a(e). Thus, de^{e) is 
given by: 



1 if /ce (e) < 0, 

(9e</? (e) = { Pc exp (/3ce) if k^ (e) = 0, 
Pi (e) if k, (e) > 0. 



(68) 



According to (|68|l . f>(() is an piecewise monotonic 
function of e, and therefore, there is a bijective correspon- 
dence between e and 1^3(6). Hence, e and Lp are equiva- 
lent representations, where k^p = dipPip is negative; and 
therefore, the canonical-like description which uses the 
representation TZ^ — ($, N) is equivalent to the micro- 
canonical description in the thermodynamic limit. It was 
shown in this way that it is always possible to choose a 
representation where the ensemble equivalence is guaran- 
teed. 
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D. An example: the Antonov problem. 

As example, let us consider the classical Antonov prob- 
lem the self-gravitating gas interacting by means 
of the Newtonian potential, which is enclosed with a 
spherical rigid container of radio R. Considering that 
G = M = R = h = 1, the microcanonical description of 
this model system is derived from the Planck density: 



Pcif3,p)^-p\ni2TTf5)+p\np-p, (69) 

by using the thermodynamic formalism presented in sub- 
section HTUl The internal energy density e and particles 
density p are given by: 




2p 



(70) 



where C = (27r/3)~^ exp {—p) and $ = ~(3<j), being (j) the 
Newtonian potential. The numerical study is carried out 
by solving the Poisson-Boltzmann equation: 



A* = -47rexp (*) , 



(71) 



where = $ + In (fiC) = In {f3p) , imposing the following 
conditions at the origin: 



*(0) = v,*'(o) = 0, 



(72) 



being ip = In (/3po) a-n integration parameter related with 
the central density po of the system. The canonical pa- 
rameters (3 and p are obtained through the relations: 



vI/(l)=/3 + ln(/3C),vI/'(l) = -/3, 



(73) 



while the entropy and the total energy are obtained from 
the relations and (|33|l respectively. It is easy to see 
that all the thermodynamical variables and potentials are 
obtained as functions of the parameter -0. 

Figure shows the well-known caloric curve of the 
Antonov problem. The points of superior branch ABC 
correspond to the equilibrium states, while the others 
represent unstable saddle points. No equilibrium states 
are found for energies e < —0.335, which is related 
with the existence of the gravitational collapse. The 
equilibrium states belonging to the energetic interval 
—0.335 < e < —0.198 exhibit a negative specific heat, 
and therefore, they are not accessible by using the canon- 
ical ensemble. 

According to all the exposed in subsection above, the 
Antonov problem obeys to the thermodynamic limit H30|l . 
and therefore, we are able to choose an appropriate rep- 
resentation for the microcanonical variables in order to 
guarantee that the canonical-like ensemble associated 



FIG. 1: Caloric curve for the Antonov problem. The states 
beloging to the interval AB exhibit a negative specific heat, 
and therefore, they are not accessible by using the canonical 
ensemble. 



-1 ' 1 ' 1 ' T" 



—\ ' r- 



_i 1 I 1 I 1 i_ 



_I 1 L_ 



-0.4 -0.2 0.0 0.2 0.4 0.6 0.8 

energy 

FIG. 2: The curve i/j — e. There is a univocal correspon- 
dence between these variables for all equilibrium states {ABC 
branch) . 



with this new representation be equivalent to the mi- 
crocanonical one. 

Figure 101 shows the dependence between the thermo- 
dynamical variables e and ip. Disregarding the unstable 
branch, it is easy to see that there is a univocal correspon- 
dence between these variables for all equilibrium states. 
We can select an appropriate representation tp — ip (e) for 
the microcanonical description by rephrasing ijj, or more 
exactly, x = 6xp (-0) as a canonical-like parameter. The 
map e ^ if (e) can be established by using the transfor- 



9 



mation rule (|57|l for the canonical parameters: 



X = exp 



9e 



d(p jtp) 



(3 (ip) exp (-V') 



de (V-) 



(74) 

This map is shown at figure ||3J)- The bijective character 
of this representation change assures the reparametriza- 
tion invariance of Physics in the microcanonical ensem- 
ble. Figures Q and © show the "caloric curve" x~ f 
and the entropy function s — in the if representation 
for the microcanonical variable e. The first shows the 
equivalence between the canonical-like ensemble with the 
microcanonical ones, while the second one confirms the 
concavity of the entropy function in this representation. 




FIG. 3: The map: : e ip (e) . The bijective character of 
this transformation ensures the invariance of the microcanon- 
ical description. 



Thus, the physical observables obtained from the 
canonical-like probabilistic distribution function: 



exp 



-Nx ■ ^ [e /N^ 



(75) 

are equal to those obtained by using the microcanonical 
description in the thermodynamic limit 13U|) . The numer- 
ical solution of this model exposed above was obtained by 
considering fixed the parameter i/j = Inx, and this sup- 
position allows us to determine all the thermodynamical 
observables of this model system. Therefore, we are able 
to affirm that this numerical solution was obtained by 
using the x- canonical-like description. 

It is remarkable that this generalized ensemble allows 
us to access to all equilibrium states of the model, even in 
the cases in which these last ones exhibit a negative spe- 
cific heat. This is an important difference regarding to 
the usual canonical ensemble. Although the similar ap- 
pearance between these two descriptions, the simple rep- 




FIG. 4: The "caloric curve" x ~ 'fi- This figure reveals the 
equivalence between the canonical-like ensemble with the mi- 
crocanonical ones in the (^-representation. 




FIG. 5: The entropy function in the 99-representation. It can 
be note the concavity of this thermodynamical potential. 



resentation change of the microcanonical variables (com- 
patible with the scaling laws of the thermodynamical 
variables and potentials in the thermodynamic limit) al- 
lows the extension of a canonical-like description to those 
cases in which the usual Gibbs' canonical one fails. 



IV. CONCLUDING REMARKS 

We have reconsidered some aspects about the micro- 
canonical mean field description for astrophysical Hamil- 
tonian systems. Particularly, we derived a formalism 
which takes into account those cases is which the par- 
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tides also interact by means of short-range forces. This 
general framework allows us to perform the analysis of 
the necessary conditions for the existence of a power law 
self-similarity behavior of the microcanonical variables of 
the form (|22() , which is an extension of extensivity of the 
traditional systems. This study revealed that this prop- 
erty could be only satisfied by a reduced set of models 
where the nature of the short-range interactions leads to 
the relation p — between the pressure p and the inter- 
nal energy per volume e, being 7 certain constant. 

As already shown, the consideration of the 
reparametrization invariance of the microcanonical 
ensemble allows us to extend the applicability of a 
canonical-like description to those situations where 
the standard Thermodynamics based on the canonical 
ensemble is not able to describe, for instance, when 
the systems exhibit a negative specific heat. This 
possibility was illustrated in this paper by reconsidering 
the classical isothermal model of Antonov 5] . 

The physical significance of such generalized canonical- 



like ensemble relies on the reparametrization invariance 
of Physics in the microcanonical ensemble as well as the 
equivalence of these descriptions by the consideration of 
an appropriate thermodynamic limit. It is important to 
remak that such thermodynamic limit is derived from 
the scaling properties of the thermodynamical variables 
and potentials of a system with a large number of de- 
grees of freedom, and therefore, the same one can not be 
arbitrarily introduced. 

An interesting application field of this geometrical 
ideas could be in the Montecarlo methods based on the 
canonical weight exp {—(3E). It is well-known that these 
numerical methods are more simple than those based on 
the microcanonical ensemble, but the first ones diverge 
when the ensemble equivalence is not ensured. Therefore, 
the consideration of a generalized canonical-like weight 
of the form exp [-Nx ■ V {E/N'^^)] (where Ag is the 
scaling exponent of energy in the thermodynamic limit 
N 00), could guarantee the convergence of such meth- 
ods for all the energetic range. 
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